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In neutron capture for magic-shell nuclei the direct reaction mechanism can be important and 
may even dominate. As an example we investigated the reaction 48 Ca(n,7) 49 Ca for projectile energies 
below 250 keV in a direct capture model using the folding procedure for optical and bound state 
potentials. The obtained theoretical cross sections are in agreement with the experimental data 
showing the dominance of the direct reaction mechanism in this case. The above method was also 
^ ■ used to calculate the cross section for 50 Ca(n,7) 51 Ca. 
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I. INTRODUCTION 

> . 

Attempts to understand low- and medium-density (n n < 6- 10 16 n/cm 3 ) neutron-capture nucleosynthesis processes 
I require temperature-dependent Maxwellian-averaged capture cross sections. For radioactive isotopes, where these 
cross sections cannot be measured they are usually calculated within global Hauser-Feshbach models. Whereas this 
' statistical approach appears to be applicable for nuclides in between neutron shells exhibiting sufficiently high level 
, densities, it certainly will fail for neutron-magic isotopes with only a few widely spaced resonances. For such cases 
OS ' apart from compound-nucleus (CN) capture, direct (DI) reaction processes are important. 

As an example we investigated the reaction 48 Ca(n,7) 49 Ca to test our method, where, in principle, all necessary 
information is already available from various experiments We also calculated the reaction 50 Ca(n,7) 51 Ca, which 

may be of importance in connection with an interpretation of the 50 Ti overabundance observed in FUN inclusions of 
the Allende meteorite (see e.g. JlO| , pT| ). 

In the CN mechanism the projectile merges with the target nucleus and excites many degrees of freedom of the 
CN. The excitation proceeds via a multistep process and has a reaction time typically of the order 1CP 14 s to 10~ 20 s. 
After this time the CN decays into various exit channels. In the DI process the projectile excites only a few degrees 
of freedom (e.g. single-particle or collective). The excitation proceeds in one single step and has a characteristic time 
scale of 10 -21 s to 10~ 22 s. This corresponds to the time it takes the projectile to pass through the target nucleus; 
this time is much shorter than the reaction time of the CN processes. 

For thermonuclear (in the keV region) and thermal (in the meV region) projectile energies the competition between 
the different reaction mechanisms is quite complicated. Normally the CN formation prevails below projectile energies 
of approximately 10 to 20 MeV. However, for light nuclei and magic nuclei the CN formation is often suppressed 
because there are no CN levels that can be populated. In this case the DI reaction mechanism can dominate the 
nuclear reaction. 



II. NUCLEAR-STRUCTURE INFORMATION 



The nuclear structure of the compound nucleus 49 Ca is of special interest because of the expected simple single- 
particle character of the low-lying states formed by coupling 2p 3 / 2 , an d 2pi/ 2 neutrons to the J™ = + ground 
state of 48 Ca. At higher excitation energies many bound levels can be described by 2particle-lhole configurations ||. 
Information on neutron-unbound states in 49 Ca (S n = 5.142 MeV) comes from three different experimental sources: 
neutron-capture and transmission measurements 1^,^-^, (d,p)-reaction work and high-resolution spectroscopy 
of /3-delayed neutron (/3dn) decay @,§ of J* = 3/2+ in 49 K. As is discussed in Ref. p|, the 49 K(/3-) 49 Ca(n) 48 Ca 
decay can be considered as inverse process to s- and d-wave neutron capture on 48 Ca. Companion theoretical studies 
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to Ref. Q by Divadeenam et al. |L3| (2p-lh doorway model), to Refs. by Dobado and Poves (l4| (complete 
sd-fp shell model) , as well as own QRPA shell-model calculations (Gamow-Teller strength), using the code of Moller 
and Randrup Q agree quite well with the experimental observations. The most remarkable result is the absence 
of J 71 ' = l/2 + states (s-wave resonances) up to at least S n + 800 keV, which can be understood in terms of the 
specific quasiparticlc structure of 49 Ca. The second observation is that the d-wave strength is considerably larger 
than the p-wave strength. This is due to the existence of a d-wave giant resonance for A = 48. Finally, no strong 
neutron-resonances were identified below 158 keV. The two very small resonances at 20keV and 107 keV only found 
in neutron capture , but not in /3-delayed neutron decay ||J| and not in the (d,p)-reaction |^,^| , are probably also 
d-wave resonances with T n < T 7 . 

With this nuclear-structure information on 49 Ca average continuum (HF) and resonance (Breit-Wigner, BW) 
neutron-capture rates were derived using the code SMOKER |jl2|,^6| ■ When comparing these values with the measured 
30keV Mawellian-averaged capture cross section of 48 Ca (See Refs. ||J^]), one can draw the following conclusions: 
The statistical-model (HF) prediction agrees with the measured cross section at 30 keV of about 1 mb; this result 
must, however, be regarded as completely fortuitous. The contribution of the resonances, i.e. the BW cross section, 
is only about 5% of the total neutron-capture rate. Hence, in the case of the doubly-magic nucleus 48 Ca, 95% must 
be due to direct reaction processes. 

In contrast to Ca(n,7), the nuclear structure information on the radioactive target 50 Ca and the compound nucleus 
51 Ca are scarce. Some bound levels without spin assignment are known from reaction work and 51 K /3-decay [ p"7| , 
and the singles spectrum of /3dn's has been measured (see Fig. [j], taken from Ref. the continuum underlying the 
peaks in the spectrum is due to the response of the spectrometer used). Due to Gamow-Teller selection rules, the 
peaks seen in the spectrum correspond to J"" = 1/2+ and 3/2 + neutron-unbound states in 51 Ca. According to the 
inverse relationship between /3dn-decay and neutron capture Q , these states represent s- and d-wave resonances in 
50 Ca(n,7) 51 Ca. Similar to the compound nucleus 49 Ca, the most remarkable result is the absence of J T = 1/2+ states 
up to 0.85 MeV beyond the neutron separation energy (S n — 4.4 MeV iQ). This results in a resonance capture cross 
section of (c CN ) ~ 8.5 • 10~ 15 mb when using the BW-formalism of the SMOKER codejjl|,0. If this were the total 
capture rate for 50 Ca, in any astrophysical s- and n/3-type neutron-capture process j^jioj the build-up of A > 50 Ca 
isotopes would be strongly hindered by successful competition of /3-decay. With this, a strong overabundance of 50 Ti 
— as observed in certain meteoritic inclusions (for discussion see Refs. (100) — would result. 

However, also in the case of 50 Ca neutron capture, it is not unlikely that (c ) represents only a small fraction of 
the total cross section which may be dominated by the DI reaction rate. Shell-model considerations (using the QRPA 
code of Moller et al. [|l5|,^9)) support this possibility. When assuming the compound nucleus 51 Ca to be spherical, 
its ground state has the ^3/2 configuration (J w = 3/2~). The two lowest excited levels are predicted to be the 
vpi/2 and vis/2 shell-model states at about 1.177MeV and 1.209 MeV, respectively. From the + ground state of 
50 Ca indeed strong p-wave radiative capture can be expected to the J T = 3/2~ and l/2~ levels in 51 Ca. Another 
interesting result of our QRPA calculations is the prediction of the lowest-energy J 71 — l/2 + and 3/2 + states of vgg/2 
shell-model origin around 5.50 MeV in 51 Ca. This is in good agreement with our experimental result from the 51 K 
/3dn-spectrum (Fig. [j]). After correction for recoil, we obtain for the two lowest neutron-unbound = 1/2+ , 3/2+ 
states El = (S n + E nl ) ~ 5.31 MeV and £ 2 * = (S n + E n2 ) ^ 5.37MeV. The neutron peaks at 1.18 MeV, 1.46 MeV, 
2.20 MeV and 2.48 MeV identified in our /3dn-spectrum correspond to transitions to the 2+-state in 50 Ca and, hence, 
are not important for the present resonance neutron-capture considerations. 

In the following, we will present our direct neutron-capture calculations using as input parameters the experimental 
and theoretical information given above. 



III. DI-CALCULATIONS 



Potential models have often been used to describe direct reactions at thermonuclear and thermal projectile en- 
ergies (Ref. [^(J and references therein). They are based on the description of the dynamics of the reaction by a 
Schrodinger equation with optical potentials in the entrance and exit channel. Such models are the Distorted Wave 
Born Approximation (DWBA) for transfer reactions or the Direct Capture Model (DC) for radiative capture. 

The DC cross section ist given by [^l| 

^(Qi^iW^Jshi £ I T MA M a M B ,* I 2 • (i) 



Ma M a M B , a 
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The quantities I a, Ib and Sa (Ma, Mb and M a ) are the spins (magnetic quantum numbers) of the target nucleus 
A, residual nucleus B and projectile a, respectively. The reduced mass in the entrance channel is given by /z. The 
polarisation a of the electromagnetic radiation can be ±1. The wave number in the entrance channel and for the 
emitted radiation is given by k a and kj, respectively. 

The multipole expansion of the transition matrices Tm a M a m b ,<j including electric dipole (El) and quadrupole (E2) 
transitions as well as magnetic dipole (Ml) transitions is given by 



Tm a M a Mb, a — ^M A M a M b ,ct dg cr(fi) + ^Ma M a Mb, tr^Sai^) + ^Ma M a Mb, tr^Strffi)- 



(2) 



The rotation matrices depend on the angle between k a and k 7 which is denoted by where <5 = Ma + M a — Mb- 
Defining 



m a m A J 



J_ 2 f Za 

12 \ml m 2 , 1 ' 



(3) 



we can write for the transition matrices for the electric dipole (EC = El) or quadrupole (EC = E2) transition 

Tff A M a M B , a = Y. lla ( la0Sa Ma I ^ M «) Mb - MA lA Ma I lB Mfi ) 
1*3* 

x (CS Jb M B - M A | j a M a )C(EC) l a l b j b 



x (l b 0C0 | l a O)W{Cl b j a S a ;l a j b )I^ hlB ,i 



(4) 



In the above expressions Z a , Za and m a , tha are the charge and mass numbers of the projectile a and target nucleus 
A, respectively. The quantum numbers for the channel spin in the entrance channel and for the transferred angular 
momentum are denoted by j a and j b , respectively. Here and in the following the abbreviation I stands for v2l + 1 
(the same applies to other quantum numbers). 
For magnetic dipole transitions we obtain 



rpMl _ \ ^ 
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(l a 0S a M a | j a M a ) (j b M B - M A I a M a \ I B M E 
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< A* ( — T + — \ ) hjb \/la (la + l)W(l l a j a S a \ l a jb) 

1 m A mi 



+ 2 Ma S a j b y/Sa (S a + l)W (I S a ja h\ S a jb) 



(la 0S a M a \ j a M a ) (j a M a I A M B ~ M a \ I B M B ) 



x (I a Mb — M a IS \ Ia Ma) 



Ia + 1 
Ia 
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2 m p c 2 
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where W is the Racah coefficient, the fii are the magnetic moments and m p is the mass of the proton. 
The overlap integrals in Eqs. ([|,||) are given by 



ffitlB-Jai *= drU lb]blB (r)0^(r) XlaJa (r) 



(6) 



for the electric dipole (EC = El) or quadrupole (EC = E2) transition, and by 
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rMl 

*h jb I. 



,;l aja = J drU lbJblB (r)O ml (r) Xla]a (r) (7) 



for the magnetic dipole transition (ML = Ml). The radial part of the bound state wave function in the exit channel 
and the scattering wave function in the entrance channel is given by Ui b j b i B (r) and \i a j a ( r ); respectively. The radial 
parts of the electromagnetic multipole operators are 

M1 (r) = ±- [sin p + p cos p], (8) 
2p 

O e1 (r) = 4r \{p 2 - 2) sin p + 2 p cos pi r, (9) 
pi 

G E2 (r) = i| [(5/j 5 - 12) sin p + (12 - p 2 ) p cos pi r 2 . (10) 

p J 

In the long wavelength approximation — applicable in our case, since p = fc 7 r <C 1 — these quantities reduce to 

M1 (r) ~ 1, (11) 
E1 (r)~r, (12) 
E2 (r) ~ r 2 . (13) 

The most important ingredients in the potential models are the wave functions for the scattering and bound states 
in the entrance and exit channels. In calculations performed by our group the potentials are determined by using the 
folding procedure p0[ . In this approach the number of open parameters is reduced considerably compared to more 
phenomenological potentials (e.g. Saxon- Woods potentials). The nuclear densities are derived from nuclear charge 
distributions [p3[ and folded with an energy and density dependent nucleon-nucleon (NN) interaction v g p4| : 



V(R) =\V F (R) = A J J p a (ri)pA(r2)«cff(S,p a ,p j4 ,s)df 1 dr 2 . (14) 

The variable s in the NN interaction term is given by 

s = \R + f 2 -f 1 \ (15) 

with R being the separation of the centers of mass of the two colliding nuclei. The normalization factor A accounts 
empirically for Pauli repulsion effects and dispersive parts in the potential V(R) that are not included in the folding 
potential Vp(R). This parameter can be adjusted to elastic scattering data and/or bound and resonant state energies 
of nuclear cluster states and at the same time ensures the correct behaviour of the wave functions in the nuclear 
exterior. 

For the calculation of the DI capture cross section (Eqs. (0-0)) we used the code TEDCA psfl . The folding 
potentials (Eq. (|l4|)) were determined with the help of the code VOLD f p6[ . 

The reaction 4 °Ca(n,7) 49 Ca was calculated for projectile energies below 250 keV. For the potential 48 Ca + n in the 
entrance and exit channel we used 



V(R)=\V F (R) (16) 

where Vp is the folding potential of Eq. ( |l4| ) . The strength of the potential Aj in the entrance channel was determined 
using neutron scattering data obtained from S. The data was fitted to pure s-wave scattering with a = 0.019b 
resulting in a A; of 0.93357 and giving a volume integral of 436.9 MeVfm 3 . 

For the optical potential in the entrance channel we neglected the imaginary part potential because the flux into 
the other channels is very small. One may wonder, if the capture process would not itself produce a large damping, 
because at low energies the capture cross section is larger than the elastic cross section. For instance, at thermal 
energies the neutron-capture cross section of 48 Ca is about lb, whereas the elastic neutron cross section is about 
0.02 b. However, this can be explained by the different phase-space factors for the capture and elastic cross sections. 
For the capture cross section near threshold one can write a c « (in R/k a )(— Im/ /|/o| 2 ), where /o is the logarithmic 
derivative at a radial distance R. The elastic cross section is about constant for projectiles near threshold, whereas 
the capture cross section has the well known 1/u-behaviour. From the thermal value of the capture cross section, 
using R = 4.5 fm, — Im/o/|/o| 2 ~ 6 • 10~ 5 is obtained. Since this expression can be regarded as a measure of damping 
in the entrance channel, we see that the damping effects are rather small. Therefore, at low energies the large value 
of the capture cross section (1/f-behaviour) compared to the elastic cross section (only weakly energy-dependent) is 
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due to the phase-space factors and does not result in additional damping effects. 

For the exit channels in 49 Ca we fitted Af to reproduce the experimental neutron separation energies in 49 Ca. We 
obtained A f = 0.9648, 0.8863, 0.8417, 1.6745, 1.6014 and 1.5885 for the 2p 3/2 , 2p 1/2 , lf 5/2 , 2f 5/2 , 3p 3/2 and 3p 1/2 
state in 49 Ca, respectively. Spins, excitation energies, Q-values and spectroscopic factors for the transitions to the 
ground and excited states are listed in Table | jo) . 



IV. RESULTS AND DISCUSSION 

The theoretical cross section c th is obtained from the direct capture cross section a DC in Eq. (|l]) as a sum over 
each final state i by 

a th = £c^af c (17) 

i 

In our case the isospin Clebsch-Gordan coefficients Ci arc 1. The spectroscopic factors Si describe the overlap 
between the antisymmetrized wave functions of 48 Ca + n and the final state in 49 Ca. We obtained the values for the 
spectroscopic factors given in Table | from a recent analysis of 48 Ca(d,p) 49 Ca ||. 

The most important contributions to the direct capture cross section are given by the transitions to the ground and 
first excited state in 49 Ca with a Q-value of 5.142MeV and 3.121 MeV, respectively (Table For the transition to 
the ground state (g.s.) we obtained the Maxwellian-averaged direct capture cross section at 30keV 

((Tg h s .) = 0.791 mb, 

and for the transition to the first excited state at 2.021 MeV 

(cr^ c ) = 0.222 mb. 



exc. 



With the small contributions of the 3p 3 / 2 state at 4.069 MeV and the 3p!/ 2 state at 4.261 MeV we obtain (al 
1.04 mb for the Maxwellian-averaged direct capture cross section at 30keV. The cross section for the transitions to 
the other high-spin states in the final nucleus can be neglected (Table [Oj) . This is due to the fact that only for final 
p-states an El-transition with s-wave neutrons in the entrance channel is possible. 

Fig. |^ shows the contributions to the direct capture at a projectile energy of 30keV as a function of the radial 
distance from the target nucleus given by the integrand of Eq. (^J). It is interesting to note that the important 
contributions to the investigated capture reaction come from the nuclear surface (4.5-5 fm) as well as from the nuclear 
exterior. 

In Table III the theoretical DC cross sections cr th obtained in this work are compared to the experimental data 
<j cxp . As can be seen from Table III and Fig. || the experimental data can be reproduced excellently for the thermal 
as well as for the thermonuclear energy region by our DC-calculations. In Fig. [| the theoretical direct capture cross 
section (solid line) is shown together with the experimental data of Cranston and White jl| , Kappeler et al. || and 
Carlton et al. Q| . In this plot we find the well known 1/v behaviour of the experimental as well as the theoretical 
cross section ranging from the meV to the MeV region. 

We also calculated the direct capture cross section for the reaction 50 Ca(n,7) 51 Ca. To obtain the 50 Ca + n-potential 
we used the density distribution obtained from QRPA calculations fl5||. The Q-value for the p 3 / 2 transition to the 
ground state is given by 4.4 MeV ]l8| . The resonant states J n — l/2 + and 3/2 + of vg§/2 shell-model origin around 
5.50 MeV mentioned in the introduction are of 2p-lh type and therefore cannot be reached by direct capture. They 
have not been included in our calculation. The spin and the excitation energy of the first excited state is l/2~ at 
1.1768 MeV. The above values and the spectroscopic amplitudes were obtained from QRPA shell model calculations 
using a folded- Yukawa single-particle and a Lipkin-Nogami pairing model [Esi. These input parameters (see Table 



IV|) yield 

(a DC ) = 0.7 mb 

for the Maxwellian-averaged cross section of 50 Ca(n,7) 51 Ca at 30keV. The lower value of this cross section compared 
to the cross section of 1.04 mb for 48 Ca(n,7) 49 Ca is mainly due to the lower Q-value. 

These results are important for astrophysical r- and a-process calculations. Previously, the rates were calculated 
in the statistical model which yields a vanishing cross section for 50 Ca(n,7) 51 Ca. Therefore, the turning point in the 
neutron-capture path was at 48 Ca because of the very small cross sections of the heavier isotopes. With the above 
considerably higher cross section (comparable to the one of 48 Ca(n,7) 49 Ca) the path can continue beyond 48 Ca or even 
beyond 50 Ca. Recent network calculations prove that possibility at least for low entropies of the Type II supernova 
hot entropy bubble |27|]. This underlines the importance of direct capture in astrophysical environments. 
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V. SUMMARY 



We used the reaction 48 Ca(n,7) 49 Ca as a test case for calculating the DC cross section, since for this reaction all 
the relevant information is available from various experiments. The experimental cross section could be reproduced 
assuming a direct mechanism and using the potential model and the folding procedure. 

As a more general conclusion, we have confirmed for the 48 Ca region that the applicability of the statistical 
assumptions in the commonly used HF calculations to derive neutron-capture cross sections breaks down near magic 
shells. Rather, the cross section is dominated by one or a few resonances, or — as in our cases — by direct radiative 
capture to bound final states in the absence of low-lying CN resonances. 

As we have shown, some experimental information on neutron-capture resonances — even far from stability - 
can be obtained from the decay model of /3-delayed neutron emission. With regard to spectroscopic factors of bound 
states, they may be obtained from (d,p)-reactions in inverse kinematics using cooled radioactive beams. 
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FIG. 1. Beta-delayed neutron spectrum of J n = l/2 + 51 K decay. 

FIG. 2. Contributions to the direct-capture cross section at 30keV projectile energy as a function of the radial distance from 
the target nucleus. 

FIG. 3. Comparison of the calculated direct-capture cross section (solid line) with the experimental data of 48 Ca(n,7) 49 Ca 
for projectile energies from the meV to the MeV region [1],[6],[7]. 



TABLE I. Spins, excitation energies, Q-values and spectroscopic factors of levels in 49 Ca [9]. 


J" 


E cx [MeV] 


Q [MeV] 


S 


2P3/2 


0.000 


5.142 


0.84 


2pi/ 2 


2.021 


3.121 


0.91 


lf S /2 


3.586 


1.556 


0.11 


2f5/2 


3.993 


1.149 


0.84 


3P3/2 


4.069 


1.073 


0.13 


3pi/2 


4.261 


0.881 


0.12 




TABLE II. Maxwellian-averaged direct-capture 


cross sections for neutron capture in 48 Ca. 




,F 


(a (25 keV)) [mb] 




(30keV)) [mb] 


2P3/2 


0.872 




0.791 


2P1/2 


0.244 




0.222 


%2 


0.000 




0.000 


2f5/2 


0.000 




0.000 


3P3/2 


0.022 




0.020 


3P1/2 


0.007 




0.006 


Sum 


1.145 




1.039 


TABLE III. Theoretical and experimental cross section for the capture reaction 48 Ca(n,7) 49 Ca. 


ENERGY a 


REFERENCE OF 


EXPERIMENTAL DIRECT CAPTURE 




EXPERIMENTAL DATA 


CROSS SECTION CROSS SECTION 


0.0253 eV (M. 


a.) Cranston and White Jl] 


(1.09 ± 0.14) b 


1.13b 


25keV (M. a. 


.) Kappeler et al. JrJ 


(1.03 ± 0.09) mb 


1.15 mb 


30keV (M. a. 


.) Carlton et al. [ff 


(1.05 ± 0.13) mb 


1.04 mb 


97keV 


Kappeler et al. [0] 


(0.55 ± 0.09) mb 


0.58 mb 


a M. a. Maxwellian averaged. 








TABLE IV. Spins, excitation energies, Q-values and spectroscopic factors of levels in 51 Ca [15]. 




r 


E cx [MeV] 


Q [MeV] 


S 


2P3/2 


0.000 


4.400 


0.80 


2pi/ 2 


1.177 


3.223 


0.92 
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